Spontaneous P-violation in QCD in extreme conditions 

A. A. Andrianov* D. EspriuJ 
Departament d'Estructura i Constituents de la Materia and 
ICCUB Institut de Ciencies del Cosmos, 
Universitat de Barcelona, Diagonal 647, 08028 Barcelona, Spain 
V.A. Andrianov* 

V.A. Fock Department of Theoretical Physics, St. Petersburg State University, 
ul. Ulianovskaya 1, 198504, St. Petersburg, Russia 

Abstract 

We investigate the possibility of parity being spontaneously violated in QCD at finite baryon density 
and temperature. The analysis is done for an idealized homogeneous and infinite nuclear matter where the 
infiuence of density can be examined with the help of constant chemical potential. QCD is approximated 
by a generalized a model with two isomultiplets of scalars and pseudoscalars. The interaction with the 
chemical potential is introduced via the coupling to constituent quark fields as nucleons are not considered 
as point-like degrees of freedom in our approach. This mechanism of parity violation is based on interplay 
between lightest and heavier degrees of freedom and it cannot be understood in simple models retaining 
the pion and nucleon sectors solely. We argue that, in the appropriate environment (dense and hot nuclear 
matter of a few normal densities and moderate temperatures) , parity violation may be the rule rather than 
the exception and its occurrence is well compatible with the existence of stable bound state of normal nuclear 
matter. We prove that the so called 'chiral collapse' never takes place for the parameter region supporting 
spontaneous parity violation. 

Preprint UB-ECM-PF/09-09; ICCUB-09-188 



1 Introduction 

The appearance of parity (P) violation via pseudoscalar condensation for sufficiently large values of temperature and/or 
chemical potential has been attracting much interest during last decades to search it both in dense nuclear matter (in 
neutron/quark stars and heavy ion collisions at intermediate energies) and in strongly interacting quark-gluon matter 
( "quark- gluon plasma" in heavy ion collisions at very high energies). At finite baryon density it was conjectured by 
A. Migdal in [T] long ago (and revisited in [2]). One should also mention the possibility of (C)P-parity violation in 
meta-stable nuclear bubbles created in hot nuclear matter [3]. Finally P violation might conceivably accompany the 
transitions to open color phases[5 such as CFL (color-flavor locking) or SC (superconducting), but these are phases 
beyond the range of validity of our analysis. While it was argued in [5] that parity, and vector flavor symmetry could not 
undergo spontaneous symmetry breaking in a vector like theory such as QCD, the conditions under which the results of 
[5] hold (positivity of the measure) are not valid for non-zero chemical potential. 

Parity violation in QCD would lead to rather remarkable experimental signals such as the same in-medium resonance 
being able to decay into even and odd number of pions, the presence of additional Goldstone bosons (six right at the 
phase transition in the exact chiral limit, and flve throughout the broken parity phase), changes in the nuclear equation 
of state, and isospin breaking effects in the pion decay constant and substantial modiflcation of the weak decay constant 
F„i for massless charged pions, giving an enhancement of electroweak decays. 

In this work we shall explore further the interesting possibility of spontaneous parity violation employing effective 
Lagrangian techniques in the range of nuclear densities where the hadron phase persists and quark percolation does not 
occur yet. Our effective Lagrangian is a realization of the generalized linear a model, but including the two lowest lying 
resonances in each channel, those that are expected to play a role in this issue. As it will hopefully become clear later this 
is the minimal model where this interesting possibility can be realized. The use of effective Lagrangians is also crucial 
to answer the second question of interest, namely how would parity violation originating from a finite baryon density 
eventually reflect in hadronic physics. 
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This work is a continuation of the one published by two of the present authors in [6]. Here we extend the analysis 
previously done to the case of non-zero temperature showing that the parity violation phase persists in some finite domain 
in the i-i — T plane. How far up in T this domain extends is not clear yet, as eventually one is expected to enter in a 
deconfined phase for which the effective meson-quark Lagrangian is not valid. 

We have considered some of the thermodynamic properties of dense baryon matter in our approach and we have 
also addressed the issue of how our model can describe the saturation point and the formation of stable nuclear matter. 
We had not touched upon this issue in detail in [6]. We find that our description turns out to be rather accurate in 
describing nuclear matter formation avoiding the unacceptable 'chiral collapse' [Tj usually present in quark models. 

We have also examined departures from the chiral limit, i.e. allowing for non-zero quark masses. This leads to rather 
interesting results as in this case the usual pions are not exactly massless, but the new Goldstone bosons appearing at 
the transition point to the parity violating phase are. Strong interaction phenomenology becomes indeed very unfamiliar 
at that point. These results along with other considerations shall be summarized in a separate publication [5]. 

Many techniques have been used to study QCD in unusual conditions: from meson-nucleonPji2|i9i[T0] or quark-meson 
[111 112) lagrangians for low-dense nuclear matter to models of Nambu-Jona-Lasinio type |13l 114) for high-dense quark 
matter [3]. However, for reasons explained below, all hadronic models lack for one reason or another some essential 
ingredient. One should also mention the extensive lattice work, plagued with technical difficulties when fJ. ^ [15]. Let 
us finally comment that the range of intermediate nuclear densities (from 3 to 10 times the usual nuclear density) where 
we expect parity breaking to occur is of high interest as it may be reached both in compact stars [12j and heavy-ion 
collisions [16) . 



2 A generalized sigma model for QCD 

The simplest hadronic effective theory is the linear (t- model of Gell-Mann and Levy[T7], which contains a multiplet of 
the lightest isoscalar a and isotriplet pseudoscalar tt" fields. Spontaneous chiral symmetry breaking emerges due to a 
non-zero value for {a) ~ {qq) / F^. Current algebra techniques indicate that in order to relate this model to QCD one has 
to choose a real condensate for the scalar density, with its sign opposite to current quark masses, and avoid any parity 
violation due to a v.e.v. of the pseudoscalar density. The introduction of a chemical potential does not change the phase 
of the condensate and therefore does not generate any spontaneous parity violation. 

Thus too simple phenomenological models retaining only the lightest degrees of freedom are not capable to explore 
all the different phases that the presence of manifest CP violation due to the non-zero chemical potential opens. The 
need to include more resonances can be qualitatively motivated by the fact that at substantially larger densities typical 
distances between baryons are shrinking considerably and meson excitations with Compton wave-lengths much shorter 
than the pion one start playing an important role. 

The minimal model having the possibility of describing spontaneous parity breaking (SPB) contains two multiplets 
of scalar/pseudoscalar fields Hj — ajl + iirj, j — 1,2, with -kj = 7!""t" where r" are Pauli matrices. We require an 
exact SU{2)l x SU{2)fi symmetry in the chiral limit. We should think of these two chiral multiplets as representing the 
two lowest-lying radial states for a given J^'^. Adding more resonances would surely provide a better description of the 
hadronic phase of QCD, but the present model already possesses all the necessary ingredients to study SPB. Inclusion 
of higher-mass states is required at substantially larger densities. Likewise we do not need to include vector resonances 
(with one exception to be discussed below). 

The effective potential of this generalized a model 

+i\,{HlH2 + HlH,)HlH2^ + O(^) (1) 

contains 9 real constants. QCD bosonization rules imply that they are ~ Nc. The neglected terms will be suppressed at 
least by inverse powers of the chiral symmetry breaking (CSB) scale A ~ 1.2 GeV. If we assume the v.e.v. of Hj to be 
of the order of the constituent mass 0.2 0.3 GeV, it is reasonable to neglect these terms. 

We take Hi as the chiral multiplet coupling locally to the quark fields (see section [4} . Using the global invariance of 
the model we can parameterize 

Hi{x)=ai{x)f{x)^ai{x)e^p(^i^^ ■ ff2(a;) = ('^2(2;) + i7^2(x))c(a;). (2) 

The parities of a2(x) and n2 are even and odd, respectively (in the absence of SPB). The potential ([TJ reads 

2 ^ 

Voff = - ^ O-jAjfcCTfc - A22(7r2)^ + A2(^(7r2)^^ + (712)^ (^{X.i - \4)ai + \6a1a2 + 2X2a2^ 

+\iaf + \202 + (A3 -I- A4)(Jia2 -I- A5(Jia2 -|- A6(Tict|. (3) 
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The corresponding gap equations are 



2(Aiicri + Ai2(72) = 4Ai(7i+3A5a?(72 + 2(A3 + A4)cricr2 + A60-| + p2(^2(A3-A4)cri+Ae(72j (4) 

2(Ai20-i + A220"2) — AsCTi + 2(A3 + A4)(7iCr2 + 3A6CriCr2 + 4A20-2 + (^AeCTi + 4A2Cr2^ , (5) 

= A22 + (A3-A4)a? + A6Cri(72 + 2A20-^ + 2A2P^), (6) 

where the following notation has been introduced: (tti) = {'K'^)d^°', {1^2) — p5^°' . 



The above effective potential must exhibit the usual chiral symmetry breaking pattern at p = T = 0. For this to 
happen (cti) must acquire a real and positive v.e.v. to agree with current algebra considerations. Note that (tt") does 
not appear at all in the gap equations and hence its value is completely undetermined, but from ([2} we see that a non 
zero value for (tt") would affect the phase of the scalar condensate. The addition of a small mass for the quarks fixes the 
phase of the breaking without having to appeal to other arguments 

The previous set of gap equations may have several solutions for a\ and 02, but since we know that in normal 
conditions QCD does not break parity, p must vanish. Since for the potential to be well defined requires A2 > 0, a 
sufficient condition for the absence of SPB is 

(A3 — A4)cri + \(,(J\<J2 + 2X202 > A22. (7) 

On the other hand, the mass of the pseudoscalar tv2 is governed by the second variation 

V^^i^ = 2(-A22 + (A3 - Xi)af + AecTicra + 2A2cr^ + 6X2p^). (8) 

Positivity of this mass for p = implies ((TJ. The condition is therefore necessary too. 

Let us establish the necessary and sufficient conditions for CSB to take place in normal conditions. The necessary 
condition to have a minimum of Vefi for non-zero aj (for vanishing p) can be derived from the condition to get a local 
maximum (or at least a saddle point) for zero aj. This extremum is characterized by the matrix —Aij in ([3| which must 
have at least one negative eigenvalue. The sufficient condition follows from the positivity of the second variation of Via 
for a non-trivial solution of the two first equations @, © at p = 0. 

The matrix containing the second variations in the scalar sector is l/'^-* 

^Vj^i^ = — All + GAlCT? + 3A5(7lCT2 + (A3 + A4)o-|, 

= --2Ai2 + 3A5<T?-f4(A3 + A4)aio-2 + 3A60-|, 
^V}^J,^ = -A22 + (A3 + A4)a? + 3A6aio-2 + 6A20-|. (9) 

The required conditions are given by tr |t/(2)| > and Dett>'2' > 0, where V^^^ is the matrix of second variations 

restricted to the scalar fields. For positive matrices this means Vi^j^ > 0. The eigenvalues of ([9j eventually give the 
masses squared of the scalar mesons and thereby must be positive. 



3 The landscape of solutions 

The above gap equations exhibit a rather complex landscape of extrema and it is important to simplify the discussion as 
much as possible. Thus we shall consider here the case where the eigenvalues of the matrix — A^^ in Q are all negative. 
This corresponds to the origin being a local maximum which appears as the most natural condition. For the sake of 
searching for extrema the effective potential can be further simplified by making a general linear transformation on the 
Hj fields 

H, = J2 ^^"H^- (10) 

fc=l,2 

The Ljk must be real to preserve the reality of aj and ttj . This transformation has four real parameters which are enough 
to set Ajfc = A(5jfc and A5 = (for a more detailed discussion see [8]). Thus the particular choice of A12 = A5 = does 
not make the analysis of the extrema less general. Apart from A5, which disappears, the other effective coupling get 
modified by the above linear transformations; Xi ~* Xi . To avoid complicating the notation we shall omit the tilde for 
the transformed fields and constants in the subsequent formulae. 

One of the solution of the simplified equations is 02 ~ Q and (j\ — A/2Ai. This is a minimum, as it follows from (jS]) 
and ((9]), provided that A3 ± A4 > 2Ai > 0. Note that the last inequality follows from the stability of the potential. When 
(T2 = the effective potential exhibits an enlarged symmetry Z2 x Z2 as As does not contribute. 

If (72 7^ 0, by combining the first two gap equations it is possible to determine another set of solutions. They are 
given in term of the ratio x = 02! (y\ by the solutions to the cubic equation 

2Ai - (A3 + A4) - '^Xi.x + (A3 + A4 - 2A2)x2 ^Aex^ = 0. (11) 
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This equation may have one or three real roots. Taking into account that Ae must be taken negative (see below) and 
that the zero-order term is negative and the term linear in x is positive we conclude that there are either two positive 
and one negative real root or just one negative real root. 

In order to determine the nature of these extrema we need to compute the matrix of second variations. After solving 
for A using the gap equations and taking into account that (T2 7^ the condition for a local minimum reads 

4Ai(ai)2-iA6^ >0, 

4(A3 + A4)(Ticr2 + SAscTa, 
3 2 

-A6CriCT2 + 4A2(T2 > 0, 

-2Xiol - iAeaiaa > 0. (12) 

The first element is positive provided that As < 0, the last one is also positive if 4A4(7i < — Ascra which follows directly 
from the gap equations themselves and the condition for the absence of SPB at /i = ([7]). As for the third one, this 
gives another condition 

3 

A2cr2 > --Aeo-i a; > 0. (13) 

o 

Therefore, either there is only one saddle point (negative solution) or one minima and two saddle points (one negative, 
two positive solutions). 

As the effective potential is symmetric against the transformation ctj — » ^""i, there is actually a doubling of the 
solutions. Taking into account the solutions with (T2 = this makes a maximum of eight extrema; actually nine if we 
include the maximum at ctj = 0. Comparison with chiral algebra arguments fix the physical sign of Oj thus breaking 
the degeneracy. Of these extrema only two of them at most are minima. As to which of the two minima (should they 
exist simultaneously) has lower energy, this is actually depending on the specific values of the constants. From now on 
we revert to the original set of fields and coupling constants undoing the linear transformation (|10p . 



2 CTlCTl 

y{2) ^ 

^ cricr2 
2 7r2'7r2 



4 Inclusion of temperature and baryon chemical potential 

After bosonization the baryon chemical potential ^ is transmitted to the meson sector (in the leading order of chiral 
expansion) via a local quark-meson coupling. In turn, in the large Nc limit one can neglect the temperature dependence 
due to meson collisions and assume that the temperature T is induced with the help of the imaginary time Matsubara 
formalism for Green functions - Matsubara frequencies for quarks aj„ = {2n + l)7r//3 with P = 1/kT. In the real world 
with 3 colors this is of course an approximation, but nevertheless it should be sufficient to describe qualitatively the 
interplay between baryon density and temperature, and it is the one consistent with our mean field approach anyway. 

As already mentioned we take the chiral multiplet to have local couplings with the quark fields as being Hi. The set 
of coupling constants in ([T]) allows us to fix the Yukawa coupling constant to unity. Thus fj, and T are transmitted to 
the boson sector by the term 

A£ = - {qRHiQL + qLHlqn) > -qaiq, (14) 

where qL,R are assumed to be constituent quarks. We do not include baryon fields explicitly and therefore quark matter 
and nuclear matter are indistinguishable in our approach. Of course in the conditions where unconfined quark matter 
would be present the density would be such that our model is not applicable anymore. 

After integrating out the constituent quarks the full temperature and chemical potential dependence, to the leading 
orders in chiral expansion, find their way into the gap equations. Namely (U} is modified to 

2(Aiicri + Ai2cr2) = 4Aiai + 3X5afa2 + 2{\'j, + \4)cJia2 + Kai + {2{\z - \i)Gi + \a(j2^ + 27V(Ji^((Ti, /x, /3), (15) 
with M = ^l^J , and where 

Aiai,,,P) = 1 E exp(^^n.) / ^ (^^,^ ^ ^)2^ ,2 „ ^2 - = 0, , = 0] 



oTrlF /p. cosh(/jM)+exp(-/3g) 
= 2y., -'^i cosh(/3M)+cosh(/3£) ' 

where the Fermi distribution has been introduced. A originates from the one-loop contribution to V^ff 
normalized to vanish for very large ai. All the dependence on the environment is in the function ^(cri, /i, /3). 
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Let us particularize to the zero-temperature case. At T = 



2 ~2 2 2x3/2 4, M + - 



(18) 



and 



A{gi , ^i, 13 ^ oo) ^ 29{fi - ai) ^ " d£ ^/e^~~^ = - cr? - o"? In ^ ^ (19) 

Using the gap equations © and (|15[). the value of the effective potential at its minima is given by the compact 
expression 



It should be emphasized that all the above results have corrections of O (/i'^/A^, af/A'^). 



5 Thermodynamic properties of the model 

Thermodynamically the system is described by the pressure p and the energy density, e. The pressure is determined by 
the potential density difference with and without the presence of chemical potential, dp = —dV, 

= 14ff('^j) - Vcs(o-j{fi),p{fi),f?j, (21) 

where the dependence of aj and p(/i) on fi has been shown explicitly and aj = (7j(0). The energy density is related to 
the pressure by 

e = -p + NcPQB- (22) 

The chemical potential is defined as 

^ege = Nc^i, (23) 

with the entropy and volume held fixed. The factor Nc is introduced to relate the quark and baryon chemical potentials. 
Since e is independent of p, 

d^p^N^QB. (24) 
Thus the relation between baryon density, Fermi momenta and the chemical potential is for quark matter 

1 c, ^/ 3 f 2 , \2n3/2 /r,t-\ 

This set of identities provides a functional relation between qb and ^. The pressure can also be written as 

P = .la..(^), (26) 

showing that the energy per baryon has an extremum when the pressure vanishes. Since the pressure is an increasing 
function of the density as we have seen (and it obviously vanishes at zero density), and infinite nuclear matter is stable 
(thus implying zero pressure too) the phase diagram in the p, qb plane must necessarily exhibit a discontinuity. This 
would correspond to the first order transition associated to the formation of nuclear matter at some critical value /x*. 
This implies that several solutions of the gap equations must coexist around p* . 

Our model consisting of two scalar isomultiplets is somewhat too simple in one respect. The stabilization of nuclear 
matter requires not only attractive scalar forces (scalars) but also repulsive ones (vector- mediated) [18]. Conventionally, 
the latter ones are associated to the interactions mediated by the iso-singlet vector uj meson. Let us supplement our 
action with the free uj meson lagrangian and its coupling to quarks 

= — ^LU^^LU^" + ^mtiO^LJ^ — Qujqqqj^Lj'^q, (27) 

with a coupling constant g^jqq ~ 0{l/\/Nc). After bosonization of QCD, on symmetry grounds, any vector field interacts 
with scalars in the form of commutator and therefore ujfj, does not show up in the effective potential Hj fields at the lowest 
order. However in the quark sector the time component uq interplays with the chemical potential and it is of importance 
to describe the dense nuclear matter properties. Let us assign a constant v.e.v. for this component guiqq{(jJo) = uj. Then 
one needs to compute the modification of the effective potential due to the replacement jj, ^ fi + ui = fl. The variable Co, 
and accordingly p., is dynamical and it also appears quadratically in the mass term in (|27|l which reads 



AK. = 4-^(^0^)--^^, G..%.0(i-). (28) 
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The term (I28|l supplements the effective potential (I20p : Vcb{i^) = VcH{fi) + AVuj{jl, n) . Kft(A*) should henceforth be used 
in all the previous thermodynamical formulae. The replacement — > /7 makes all expectation values depend rather on 
p,, although of course the 'physical' p is the one directly related to the density via (I25|l . fi can be determined via the 
variation of V^ff and is given by 

^ = -N.gsM = - a,ipff\ (29) 



6 The saturation point and absence of chiral collapse 

A viable model of dense baryon matter must describe the phase transition to a stable bound state at the usual density of 
infinite nuclear matter qb = Qo = 0.15-^0.16 fm~^, the so called "saturation point". This phase transition is believed to 
be of first order similar to the vapor (hadron phase) becoming saturated and turning into liquid (nuclear matter), when 
droplets of dense baryon matter aggregate to form a homogeneous nuclear liquid. 

However in simple quark models of the Nambu-Jona-Lasinio type [7] this phase transition (for vanishing current 
quark masses) leads to a chirally symmetric phase with zero dynamical mass (zero v.e.v. of scalar fields): i.e. the 
so called "chiral collapse" takes place. This is unphysical as chiral symmetry is not restored (at least not fully) in 
nuclei. Furthermore in NJL-type models the typical baryon density is substantially larger that the normal one (typically 
Qb — 2.8qo)- For this reason the NJL model cannot be a reliable guide to phase transitions in dense nuclear matter. Let 
us investigate whether our proposal avoids these pitfalls. 

The saturation point where nuclear matter forms is characterized by vanishing pressure as is the case for the vacuum. 
Given that the pressure is a non-decreasing function of fi, this can only happen if two solutions coexist and one of them, 
initially corresponding to negative pressure, takes over. This is possible in our model because we can have up to two 
minima both of them corresponding to chirally broken vacua. Let us see how this can be implemented. 

If the value of the effective chemical potential corresponding to infinite nuclear matter /x* is such that p," < a'l the 
one-loop correction from the quark loop is zero for p. < fl* and VaB is constant up to that point. Consequently, p = 
throughout this phase. Let us now assume that there is another minimum that for fi = has higher energy that the 
previous one, hence negative pressure, and that this minimum is characterized by a value of ai < ctJ. Then chemical 
potential corrections will start modifying this second solution as soon as ^ > ai. The pressure can only increase with ^ 
and thus it will therefore cross the p = line for a value p — p* . At this point the second solution takes over and p* 
corresponds to the saturation point. Let us denote by al the value of cti at p* obtained from this second solution. 

Let us now prove that for the large set of coupling constants eventually leading to SPB (see below) one of the fields 
at least has a non-zero expectation value, aj 7^ in the chiral limit, and chiral collapse is impossible. Indeed, suppose 
that (Tj = 0,Vj at /i* . Then the matrix of second variations for the effective potential given in © and taken from psp 
reads 

iyjfli = - An + AA/i^ V}^ = -2Ai2, ^V^'X = -^22. (30) 

In order to induce spontaneous breaking of parity one has to choose A22 > (see below) . Then from (|30p one finds that 
for any signs of other constants An, A12 and for any value of /i the second variation matrix is never positive definite and 
one finds a saddle point or a maximum at the presumed saturation point and beyond deep in the nuclear matter phase. 
On the contrary, for A22 < one can always get a large enough /i so that det > and chiral 

collapse eventually takes place. As we have to guarantee the existence of stable nuclear matter we assume A22 > from 
now on. Thus there seems to be a direct relation between the absence of chiral collapse in the model and the eventual 
presence of a SPB phase. Simplifying the model to have just one isomultiplet (equivalent to NJL) simply corresponds 
to considering the first element of the matrix of second variations and this changes sign (implying chiral symmetry 
restoration) for /i large enough. 

The energy crossing condition p — at p* < ai, a* = aj{p*) can be written, taking into account (|20p and (|28p 

(a^A^fca^ - a;A,,a,!) = ^^p*p%(p*) + G^^^pUp') = Qb{h1 + G^Nlel{^i*), (31) 

where p* is related to the physical value of /i* by Ea. (|29[> . This relation represents the condition for the existence of 
symmetric nuclear matter. It can always be fulfilled by an appropriate choice of Gc^. 

This is still a crude model and one certainly should not expect an extended sigma-quark model to characterize normal 
nuclear matter with very good precision. Nevertheless taking this at face value we can derive a relation between ^* and 
al: pF = ^{p*)^ - (cri*)2 = 1.3 ^ 1.4 fm~^ = 260 ^ 270 MeV that corresponds to ^0 = 0.15 ^ 0.16 fm'^ Not bad. 



7 The SPB phase transition 

We shall consider from now on the solution corresponding to the most stable minima for p> p* . To simplify the notation 
we shall also ignore the difference between ^ and the effective chemical potential p. 
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The possibility of SPB is controlled by the inequality In order to approach a SPB phase transition when the 
chemical potential is increasing we have to diminish the l.h.s. of inequality (Ul and therefore we need to have 

[(A3 - X4)al + Ascricra + 2A2cr^] < 0. (32) 

This is equivalent to 

(Asm + 4A2a2) Vjfi, < (2(A3 - A4)ai + Agaa) F^'l, . (33) 

This last inequality is a necessary condition that has to be satisfied by the model it to be potentially capable of yielding 
SPB at large densities. 

Let us examine the possible existence of a region of fi where /o 7^ 0. Then 

(A3 - Xi)aj + Aeaio-a + 2A2 + p^^ = A22, (34) 

After substituting A22 from (I34[l into the second Eq.Q one finds that 

AsCTi + 4A4(7i(j2 + Aelcrl + p^) = 2Ai2, (35) 



where we have taken into account that ai 7^ 0. Together with p4|l this completely fixes the relation between v.e.v.'s of 
the scalar fields (Ti,2 throughout the SPB phase independently of ^ and p. If A2A6 7^ (|34l) and (|35|) allow us to get rid 
of the v.e.v. p and 

(2A5A2 + A6(A4 - A3))o-? + (8X2X4 - Xl^aia2 = 4A2A12 - A6A22, (36) 

whose solution for /j, > ficrit is 

. , B 2A5A2 + A6(A4-A3) „ _ A6A22 - 4A2A12 

+ A= JT^^^ > B= ^2_^^^^^ ■ (37) 

Let us now determine the critical value of the chemical potential, namely the value jicrit where p{ficrit) = 0, but 
Eqs.dHll, 135}, ((321) hold. Combining the two equations (|34]), ((35]) 



(4A2A12 - A6A22)a;^ + (2A6A12 - 4A4A22)a; + 2(A3 - A4)Ai2 - A5A22 =0, x=—. (38) 

In order for a SPB phase to exist this equation has to possess real solutions. If 4A2A12 — A6A22 ~ there is only 
one solution corresponding to a second order transition, but there may exist other critical points that fall beyond the 
accuracy of our low energy model (which becomes inappropriate for small values of ai). We stress that Eqs. (|37p and 
(|38[) contain only the constants of the potential and do not depend on temperature and chemical potential manifestly. 
Once we find Xcrit ~ a;±(Ai2, A22, A2, . . . , Ae) one can immediately calculate 



af{Ajt,Xj) = \j a2{/^jk,Xj) ^ x±af . (39) 

After substituting these values into Eq. (llSp one derives the boundary of the P-violation phase 

NA(at,p.,P) = All - 2Ai(a±)' - X^a^af - (A3 - Xi){<j^ f , (40) 



which is a positive combination. The relation (|40|l defines a P-breaking divide line in the T — p plane. From (|16p one can 
obtain that ^ > and ^ ^ 00 when T, p 00. It means that for any nontrivial solution x±, a^, J^A{cr^, p,j3)>Q 
the P-breaking phase boundary exists. If the phenomenon of P-violation is realized for zero temperature it will take 
place in a domain involving lower chemical potentials but higher temperatures. 

8 Nature of the SPB phase and physical spectrum 

Once a condensate for 7r2 appears spontaneously the vector SU{2) symmetry is broken to U(l) and two charged excited it' 
mesons are expected to possess zero masses. For simplicity let us consider zero temperature. After the second variations 
of the potential we define the following quantities (one factor of p is included in the definition of V for each derivation 
w.r.t. for convenience 



ivii = -Aii+6Aia?+3A5cria2 + (A3 + A4)cr2 + (A3-A4)p^+AA 



„ 2, \/p 
3(Ti In - 



V12 = -2Ai2 + 3A5a? 4(A3 + A4)aia2 + 3A6a2 4- A6P^ 

^V22 = -A22 + (A3 + A4)a? + SAgaio-z + 6A20-^ + 2A2p^ (41) 

Vio = 2p((A3 - A4)o-i +2A6(J2), 

V20 = 2p(A6Cri -I- 8A2Cr2), 

Voo = 2(-A22 + (A3-A4)a? + A6<Jio-2 + 2A20-^+6A2p'). (42) 
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The last element is equal to 8A2p^ in the region where SPB takes place. By taking one derivative w.r.t. /i of the gap 
equations and solving for fpO-j and di_ip we find 

a A\r r~'> 2 ^22Voo — V|o ..r rZ ^ V10V20 — V12V00 

^^,CJ^ = -mo^^^P -ol , ^^,(J2 = -mo^^^P -oi ^-^^^ , (43) 



Q AKf !■> 2 1^12 1^20—^101/22 

^^,p=-^N(J^^p? -gI , (44) 

At — + cxj the v.e.v have the following asymptotics: ai 0-2 ~ const, but these hold in a region where our effective 

theory is not reliable anymore. 

Let us compare the derivatives of the dynamic mass ai across the phase transition point. Their difference reads 

1 Merit +iO !Mcr,t-iO V DetVDetVj 'p^o 

provided that the determinants are positive (they determine the spectrum of meson masses squared) . Thus the derivative 
of the dynamic mass is discontinuous and the phase transition is of the second order. 

We notice that convexity around this minimum implies that all diagonal elements are non-negative. This gives 
positive masses for two scalar and four pseudoscalar mesons, whereas the doublet of charged of tt' mesons remains 
massless. Quantitatively the mass spectrum can be obtained only after kinetic terms are normalized. We just note 
that in the SPB phase the situation is rather peculiar: pseudoscalar states mix with scalar ones. In particular, the 
diagonalization of kinetic terms is different for neutral and charged pions because the vector isospin symmetry is broken: 
SU{2)v — > U{1). This triggers a rather exotic mechanism of isospin breaking via different decay constants. Even in the 
massless pion sector the isospin breaking SU{2)v — > U{1) occurs: neutral pions become less stable with a larger decay 
constant. We refer the reader to [B] for details. 

SPB also induces mixing of both massless and heavy neutral pions with scalars. In fact in the SPB phase parity is 
no longer a conserved quantity in strong interactions, so the distinction between scalars and pseudoscalars is immaterial. 
This is why while the global broken symmetry at the point of transition to the SPB is a vector one, the two Goldstone 
bosons are apparently pseudoscalars, but as emphasized the distinction is purely semantic once parity is broken. 



9 A particular example 

So far there is no sufficient experimental information to be able to fully determine the value of the nine low energy 
constants appearing in VeS. As we have seen in section |3] we can always, by means of a field redefinition, eliminate A5 
and diagonalize Aij. While using this freedom is quite convenient to discuss the landscape of extrema of the effective 
potential, this transformation does have some effects however in the coupling (|14[l to quarks qaiq g(ai +ja2)q, which 
complicates the discussion when introducing the chemical potential. 

However, just to see that the emergence of the SPB phase is quite plausible, let us assume that the low-energy model 
of QCD (m indeed has A5 = and a diagonal matrix Aij whose eigenvalues need not be degenerate or, equivalently, we 
take 7 = for the above vertex. Then following the procedure outlined in section |3] we can collect the set of inequalities 
ensuring the existence of the two minima necessary for the emergence of the saturation point 

Ai,2 > 0, Ao < 0, All > 0, A22 > 0, (A3 ± A4)Aii > 2A1A22, (A3 + A4)A22 > 2A2A11, xi > max(-^, -|^) (46) 

Ae o A2 

If in addition we want to have SPB we have to require that A3, A4 are positive (see [6]). The sign of ai is fixed by current 
algebra considerations thus removing the degeneracy of solutions. 

Both sets of solutions, namely these with a"2 = or these with (T2 7^ may provide the lowest minimum at /i = 0. 
This turns out to be controlled by the sign of the combination 8A2 A4 — Ag . The latter follows from the combined analysis 
of the vacuum solutions and the solutions at the SPB phase transition (see Sec. [T}. It is clear that the conditions leave 
a lot of room for the simultaneous occurrence of both phenomena. In [6] a rather rough phenomenological fit was done 
leading e.g. to Ai ~ 0.15, A3 ~ 4, An ~ 0.03 GeV~^, A22 ~ 0.1 GeV~'^ which would trigger SPB at about three normal 
nuclear densities. Nuclear matter is formed at /x« ~ 303 MeV stabilized by a meson condensate with Gui ~ (10 15) 
GeV"'^ in satisfactory agreement with other model estimations [19j . 



10 Conclusions 

Let us summarize here our main findings. Parity violation seems to be quite a realistic possibility in nuclear matter 
at moderate densities. We have arrived at this conclusion by using an effective Lagrangian for low-energy QCD that 
retains the two lowest lying states in the scalar and pseudoscalar sectors. We include a chemical potential for the quarks 
that corresponds to a finite density of baryons, implement the bound state of normal nuclear matter and investigate 
the pattern of symmetry violation in its presence. We have found the necessary and sufficient conditions for a phase 
where parity is spontaneously broken to exist. In general this phase is bound and it extends across a range of chemical 
potentials that correspond to nuclear densities where more exotic phenomena such as CFL or CS may occur. 
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Salient characteristics of this phase would be the spontaneous violation of the vector isospin symmetry SU{2)v down 
to U{1) and the generation two additional massless charged pseudoscalar mesons. We also find a strong mixing between 
scalar and pseudoscalar states that translate spontaneous parity violation into meson decays. The mass eigenstates will 
decay both in odd and even number of pions simultaneously. Isospin violation can also be visible in decay constants. 

We think that our conclusions are drawn in a region of parameters where effective Lagrangian techniques are applicable 
and, while obviously we cannot claim high accuracy in our predictions, we are confident that the existence of this novel 
phase is not an spurious consequence of our approach but a rather robust prediction. It would surely be interesting to 
investigate how this new phenomenon could possibly infiuence the equation of state of neutron stars . 

One could hope that lattice methods [15] may shed some light on this issue and confirm or falsify the existence of 
this interesting phase in dense nuclear matter. Work along this lines is proceeding [20] . 
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